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Abstract 

A quasiparticle model of the quark-gluon plasma is compared with lattice QCD data for purely 
imaginary chemical potential. Net quark number density, susceptibility as well as the deconfinement 
border line in the phase diagram of strongly interacting matter are investigated. In addition, the 
impact of baryo-chemical potential dependent quasiparticle masses is discussed. This accomplishes 
a direct test of the model for non-zero baryon density. The found results are compared with 
lattice QCD data for real chemical potential by means of analytic continuation and with a different 
(independent) set of lattice QCD data at zero chemical potential. 
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I. INTRODUCTION 



Strongly interacting matter, as described by QCD, exhibits an astonishingly rich phase 
structure. In the region of not too large baryon densities, the deconfinement transition from 
hadronic matter to a plasma built of quark-gluon constituents is the most prominent feature. 
It is signalled by a rapid change in the expectation value of the Polyakov loop and the chiral 
condensate where one assigns a pseudo-critical temperature T^. to this transition (cf. reviews, 
e. g. [1]). At higher temperatures, T > STc, further structural changes are conjectured 
[2]. For non-zero quark chemical potential corresponding to a finite net baryon density, 
many researchers argue on the change of the deconfinement border line, representing an 
analytic crossover, into a first-order transition curve. The onset of this sequence of first- 
order transitions is marked by a critical point being of second order which has attracted 
much attention recently (see [3]). The interest in this part of the phase diagram is triggered 
by the possibility to probe it under laboratory conditions in relativistic heavy-ion collisions. 

With the advance of precision data from ultrarelativistic heavy-ion collisions at RHIC, 
the paradigm on the quark-gluon plasma has changed [4] : The notion of a strongly coupled 
plasma has been put forward to explain the seemingly very small viscosity to entropy ratio 
deduced from hydrodynamical fits to experimental data as in [5], and various models have 
been developed [6] to account for such a property. On the other hand, we are witnessing a 
vast progress in first principle calculations of thermodynamic properties of hot deconfined 
strongly interacting matter based directly on QCD [7-9]. While various observables such 
as pressure, energy density or numerous susceptibilities are addressed, the available lattice 
QCD data are obtained for different numerical set-ups, lattice sizes, fiavor numbers and 
quark masses as well. Particular attempts are needed to access non-zero baryon densities 
because the notorious sign problem of the fermion determinant prevents a direct applica- 
tion of methods useful for zero baryon density. Nevertheless, a few methods have been 
developed to access non-zero baryon densities. Among such methods is the calculation of 
thermodynamic quantities at purely imaginary chemical potential. Here, the sign problem 
is avoided but the results have to be analytically continued to real chemical potential. In 
this respect it is useful to have a model at our disposal which is successfully probed for both, 
real and imaginary chemical potential, in order to accomplish the translation of results from 
imaginary to real chemical potential. 
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While baryon density effects are small for heavy-ion collisions at top-RHIC energies and 
will be even smaller for LHC energies, at least in the mid-rapidity region, for CERN-SPS 
and upcoming FAIR energies they are sizeable. In this respect, a firm knowledge of ther- 
modynamic bulk properties of strongly interacting matter is highly desirable. As a step 
towards achieving this goal we are going to extend our quasiparticle model [10-13] to imag- 
inary chemical potential. Here, information is obtained for /x^ < allowing, in principle, for 
identifying n — ±iijLi. The model has been tested successfully for real chemical potential [14], 
say in describing the Taylor expansion coefficients of the pressure as a series in powers of 
fi/T. In such a way, Peshier's flow equation [11] is tested in some detail. This flow equation 
transports information about the effective coupling, G^, from the temperature axis to non- 
zero n and determines to a large extent the dependence on n and thus on baryon density. 
Another important piece of the model is the quasiparticle ansatz for dynamically generated 
effective masses of quarks oc (T^ -|- When going to purely imaginary chemical po- 

tential II — > ±i// the sign of the //^ term is flipped, as also signs in Peshier's flow equation 
are changed. Therefore, the /i dependence of the model is directly tested by considering an 
imaginary chemical potential. 

In the following, two symmetries of the QCD partition function Z[T,n) are of relevance: 
(i) Z{T,iJ,) = Z{T, -/x), and (ii) Z{T,iiii) = Z{T,i{iii + ^T)), i. e. Z{T,ii) is periodic in 
/li with period 27rT/3 [15]. Symmetry (i) makes Z an even function of /i (meaning that in 
a Taylor series expansion only even powers of /i/T, and thus also of /li/T, are encountered) 
such that we can focus on /i — +ifj,i only, while (ii) implies the Roberge- Weiss periodicity [15] 
which is anchored in the center symmetry. This periodicity is characterized by lines of first- 
order transitions {Z^ transitions) at /ij = |T(1 + 2k) for all integers k and sufficiently high 
temperature T while for smaller temperatures the behavior of thermodynamic quantities is 
analytic. The endpoint of first-order transitions, T^, is determined by the crossing of the 
Roberge- Weiss transition line with the chiral critical line which is also a first-order transition 
line for Nf — A degenerate quark flavors [16] . The Roberge- Weiss periodicity implies that in 
the T — iii/T plane all sectors between jii/T = ^k and fXi/T = ^{k + 1) are copies of the 
sector between ^i/T = and fii/T = Furthermore, the subsector between ^i/T = 7r/3 
and fii/T = 2n/3 is an reflected copy of the subsector between fii/T = and /i,/T = 7r/3 
mirrored at the first Roberge- Weiss transition line. As thermodynamic quantities behave 
non- analytically at /li/T — tt/S (Roberge- Weiss transition), an analytic continuation of 
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results obtained for imaginary chemical potential to real ji has direct access to the region 
H < Tr/STonly. 

Due to the severe approximations made when linking our phenomenological model [10, 
11, 13] to QCD [12], the Roberge- Weiss periodicity is not longer apparent. Having this in 
mind, we translate the model to imaginary chemical potential in section II. The comparison 
with lattice QCD data at imaginary chemical potential is performed in section III, where 
also the continuation to real chemical potential is presented. This allows, in addition, for a 
comparison with another and independent set of lattice QCD data obtained at // = (sec- 
tion IV). Furthermore, we investigate in detail the impact of the baryo-chemical potential 
dependence of the quasiquark and quasigluon masses (selfenergies) on the found results and 
discuss the deconfinement border line in the phase diagram of strongly interacting matter. 
Our results are summarized in section V. Appendices A and B contain Peshier's flow equa- 
tion for imaginary chemical potential and a discussion about the parametrization of the ^ 
dependence of the density. 



II. QUASIPARTICLE MODEL AT IMAGINARY CHEMICAL POTENTIAL 

The employed model is based on a two-loop $ functional approach to QCD with corre- 
sponding one-loop self-energies considered in HTL approximation in the asymptotic limit 
and the neglect of finite width effects, (anti)plasmino and longitudinal gluon contributions 
as well as Landau damping [12]. The QCD running coupling is replaced by an effective 
couphng G^{T,ijl) which is subject to Peshier's flow equation [11] resting on a thermo- 
dynamic self-consistency condition and the stationarity of the grand canonical potential 
fl = —pV = —T In Z, where p denotes the pressure and V the volume of the system. 
Straightforward replacement of fi = ifii in p{T,fi) renders the net quark number density, 
n{T,iiii) — —idp{T,iiii)/diii, related to the net baryon density ub — ^n, to 

~ n^Jo W^^/^ + 2e-^/^ cos(//,/r) + ^ ^ 

where dq — 2NcNf is the degeneracy factor of quarks for Nc — 3 colors and Nf quark 
flavors. The found result for n is purely imaginary and positive (negative) for small positive 
(negative) /li, i. e. n is an odd function in /li. 



The quark dispersion relation u!q{k) employed in Eqs. (1,2) reads 

^k' + Ml (3) 

with asymptotic mass = rrig + 2M^ using 

TV^ - 1 ip 
K = -T^iT^ - ^)G\T,ti,,) (4) 



as plasma frequency. This dispersion relation is based on a calculation of one-loop self- 
energies with finite quark masses in Feynman gauge in the asymptotic limit [17] for small 
rUg/T, where rUq may be temperature dependent as well to allow for a direct comparison 
with lattice QCD data. (A different approximation of M^o is discussed in section IIIC.) 
Eqs. (1, 2) highhght the quasiparticle character of the model: the baryon charge is carried 
by excitations with dispersion relation given by Eq. (3). The dependence of on the 
chemical potential (cf. [18]) will be discussed in section HID. 

Peshier's flow equation [11] for imaginary chemical potential reads 

& = a^^ + a,.^, (5) 

where the coefficients b, ar and a^. depending on T, //j and G^{T,iiii) are relegated to Ap- 
pendix A. Transforming Eq. (5) to a system of three coupled ordinary differential equations, 
it can be solved by the methods of characteristics knowing, for instance, G^(T, = 0). A 
convenient parameterization of G^(T, = 0) is [11] 

«^(^^^-''-°) = fl!S' 

making some contact to perturbative QCD at very large temperatures. Here, Po = ^{llNc — 
2Nf) and ^ is parametrized phenomenologically as ^ = A(T — Ts)/Tc with scale parameter 
A and Tg shifting the infrared divergence to T = -|- Tc/ X < T^. for appropriate parameters, 
while we focus here on the region T > T^. 

Results obtained for //^ < need to be analytically continued into the //^ > half-plane 
in order to achieve physical results. An effective analytic continuation requires a positive 
second derivative of Z with respect to /x, cf. [19, 20], i. e. the quark number susceptibility 
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x(T, //) = dn{T,ii)/diJi, > 0. The QPM result for x reads for imaginary chemical potential 

da 



27r2T 

dn 



oo (2e^<,/T cos{iii/T) + 4e2'-«/^ + 2e'^^/^ cos(/ii/T)) 



X 



2tt^T 



dk 



(gSo-./T + 2e^«/^cos(/Xi/T) + 1)^ 



o^^ (e2^,/T + 2e-?/^ cos(/ii/T) + 1)^ 



(7) 



t-t2 _ 

7r2 



it is purely real and symmetric under /ij —fii (cf. Appendix A). Furthermore, for small 
/li, the first term in Eq. (7) is positive and dominates the second term. At /i — 0, one finds 

d„ r e^^/^ 



/ 

Jo 



dkk' 



>0, 



(8) 



where = uiq{T, = 0). 



III. COMPARISON WITH LATTICE QCD RESULTS FOR IMAGINARY MU 
A. Baryon density and quark number susceptibility 

We confront now the above introduced quasiparticle model (QPM) with lattice QCD 
data [19, 21] at non-zero T and Hi obtained for Nf — A degenerate quark flavors with 
ruq — 0.2 T; these calculations [19, 21] are performed on a lattice with temporal and spatial 
extensions A^^- = 4 and A^. = 16. Valuable information in such simulations with imaginary 
chemical potential is obtained for fii < [22] implying a negative imaginary part of the 
net quark number density n{T,ifii) according to Eq. (2). In the following, however, we 
will consider /li > which renders the sign of n{T, i/ii) and accordingly the behavior of 
X — dn/d{ifj,i). Our model is formulated for a system infinite in space and time. Thus, 
we need a proper extrapolation of the lattice QCD data to the continuum limit (A^^ — > oo 
at fixed temperature). Different estimates for a continuum extrapolation are conceivable. 
For instance, one may select a scaling factor strictly valid only for asymptotically high 
temperatures, or one may use as scaling factor the ratio of thermodynamic quantities for a 
massless, non-interacting gas of quarks and gluons known in the continuum limit and from 
lattice QCD for finite A^- Even though such estimates for a correction factor could depend 
on T, in general, we apply the latter procedure, assuming that the continuum extrapolations 
for QCD and for the non-interacting gas of quarks and gluons are similar (cf. discussion in 
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FIG. 1: Comparison of the QPM (solid curves) for the scaled net quark number density n/T^ as 
a function of in/Tc with continuum estimates of the lattice QCD data [19, 21] for temperatures 
T = 1.1, 1.5, 2.5, 3.5 Tc (diamonds, circles, squares and triangles, respectively). The fat cross depicts 
the Roberge- Weiss critical chemical potential ^ic/T = 7r/3 for T = 1.1 where we stopped our 
calculations. 

[23] for pure SU(3) theory). In principle, however, a profound extrapolation to the continuum 
limit should be based on simulations with different size lattices as leading corrections to the 
continuum limit are of the order 0{N~'^) [23, 24]. Taking the Stefan-Boltzmann result of 
ub/T^ for A^^ = 4 [19], we find as educated guess for the needed continuum extrapolation 
factor of the net quark number density c/^^j = 0.456. This compares well with continuum 
extrapolation factors reported in [25, 26], reading 0.446 and 0.465 respectively, where similar 
actions have been used in the lattice simulations. 

In Fig. 1, we compare our model with the continuum estimate of the lattice QCD data 
[19, 21] for the scaled net quark number density as a function of ^i/T^ at constant T. Because 
n = 3 as a function of imaginary chemical potential is found to be purely imaginary, both 
in Eq. (2) and in the lattice calculations, we exhibit its imaginary part in the following. 
The parameters of the effective coupling G'^{T,fi = 0) in Eq. (6) read Tg = 0.96 Tc and 
A = 56 shifting the divergence of G'^{T,fi = 0) to approximately T = 0.98 T^, where we 
utilize Tc = 163 MeV as given in [16] for the case at hand. Note that we consider only 
temperatures T > Tc. The continuum extrapolated lattice QCD data, in particular the 
pronounced bending of n/T^ for T = 1.1 Tc, are impressively well described by the QPM 
parametrization. The drastic change in the slope for T = 1.1 Tc signals the onset of the 
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FIG. 2: Continuation of the QPM results for n/T^ exhibited in Fig. 1 to real chemical potential 
/i/Tc (solid curves) for T = 1.1, 1.5, 2.5, 3.5 Tc (from top to bottom). For comparison, we also 
show the analytically continued results (dashed curves) of the polynomial fit from [19] to n/T^^ for 
imaginary chemical potential. 

Roberge- Weiss transition at ^c/Tc = Htt/SO, where n should exhibit a discontinuity. In the 
QPM, this change in slope is driven by the dependence of the quasiparticle asymptotic mass 
Moo on chemical potential and, in particular, by the behavior of with respect to /Xj as 
dictated by Peshier's flow equation Eq. (5). We note that n/T^ exhibited as a function of 
/ij/T shows almost no dependence on T for temperatures T > 1.5 Tc. Below Tc, however, 
n/T^ displays a qualitatively different behavior being continuous and periodic as a function 
o^^^^/T [19]. 

Within the QPM, results obtained by considering purely imaginary chemical potential 
can easily be analytically continued to real /i. This is achieved by continuing the purely 
imaginary variable = i/ij to the entire complex plane and finally taking the limit Im yU — > 0. 
In this way, we recover the quasiparticle model [10-12] formulated for real /x. Within the 
analyticity domain, i. e. for /i < fidT), the analytic continuation is unique as guaranteed 
by general arguments. Keeping the QPM parameters A and T, fixed, the results of n/T^ 
for real /i/Tc are exhibited in Fig. 2 (solid curves). These results may be compared to other 
analytic continuations. For instance, in [19], a polynomial fit to n/T^ as well as its analytic 
continuation to real /i/T (dashed curves in Fig. 2) was considered. Despite the fact that this 
polynomial fit n(T, /ij, m^) = a(T, mg)/ij + 6(T, mq)fif for imaginary chemical potential, with 
analytic continuation n(T,fi,mq) = a(T,mq)fi — b(T,mq)n^, cannot account for the change 
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FIG. 3: Scaled net baryon density ub/T^ as a function of T/Tc for constant imaginary fj^s/Tc = 
3ifj,i/Tc (solid curves) and for corresponding real fis/Tc (dashed curves). Note that for all temper- 
atures, dn/dT > is fulfilled, as required from thermodynamic stability conditions. 

in slope observed for T = 1.1 Tc at large fii/Tc, its coefficients a and b are temperature 
and quark mass dependent. In contrast, the QPM parameters A and Tg are once adjusted 
to n/T^ at T = 1.1 Tc (cf. Fig. 1) and then kept fixed for all temperatures and chemical 
potentials. In addition, the behavior of analytic continuations of polynomial fits decisively 
depends on the considered order in /i^ (cf. discussion in [9, 27]). The QPM, in contrast, 
contains all orders of fif respecting the symmetry In Z{fi) = \nZ{—^). As evident from 
Fig. 2, we point out that close to ^c{T) a sensible analytic continuation is needed. 

In Fig. 3, the net baryon density hb/T^ is exhibited as a function of T/Tc for constant 
imaginary (solid curves) as well as for real baryo-chemical potential /i^ = 3/i (dashed 
curves) . Somewhat surprisingly, the results for real /i^ significantly deviate from the original 
results for imaginary chemical potential only at large /i^ and temperatures close to T^.. Note 
that in these considerations /is is restricted to \^b\ < ttT. 

Susceptibilities are quantities serving as measures of fiuctuations. The quark number 
susceptibility x (cf- Eqs. (7, 8)) is simply the derivative of the density in /Xj direction. We 
exhibit x/T"^ either at = for various temperatures (Fig. 4, left panel) or for T = 1.1 Tc 
for various values of /ij (Fig. 4, right panel). Clearly, if lattice QCD data for n{T,ifii) are 
well described by a model, the model should also describe x(T, i^i). This is indeed the case, 
see Fig. 4, where both, lattice QCD data as well as QPM results, are obtained by numerical 
differentiation of the net quark number density. The only concern that could arise is that 
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FIG. 4: Left: comparison of the QPM (solid curve) for the scaled quark number susceptibility 
x/T"^ as a function of T /T^ for /x = with the continuum estimate of the lattice QCD data in [19] 
(circles). Right: comparison of the QPM (solid curve) for x/^^ &s a function of Hi/Tc < 1 for 
T = 1.1 Tc with the continuum estimate of the lattice QCD data in [21]. 

derivatives enhance possible systematic differences between a model and the data. Fig. 4 
does not point to such a possibility. 

B. Deconfinement border line 

The solution of Peshier's flow equation Eq. (5) is accomplished by the method of charac- 
teristics. As in [11], we consider the characteristic curve emerging at T = Tc and /i = as an 
indicator of the pseudo-critical line. This transition line has been calculated in lattice simu- 
lations [16] for imaginary chemical potential. The lattice QCD data have been analyzed by 
applying polynomial fits which were analytically continued to real [16]. The results of such 
analytic continuations decisively depend on the chosen degree of the considered polynomial 
or ratios thereof, as discussed in [9, 27]. 

In Fig. 5, the phase diagram is exhibited in specific coordinate systems. Negative values 
/i^ < indicate purely imaginary baryo-chemical potential, whereas positive /i^ > indicate 
real /i^- Diamonds represent the polynomial fit [16] for both, imaginary chemical potential 
and the corresponding analytically continued results. For comparison, we depict the QPM 
characteristic (solid) curve starting at T = as solution of Peshier's flow equation Eq. (5) 
for imaginary chemical potential as well as for real /i^- The flatness of the curve in the 
exhibited /i^ interval (left panel) signals the dominance of the /z^ term in agreement with 
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FIG. 5: Phase diagram for imaginary and real baryo-chemical potential. Left: T/T^ vs. fi'^/T^. 



Right: T/Tc vs. ^b/Tc- Details are explained in the text. 

the polynomial fit findings in [16]. Note, however, that the QPM result contains all orders of 
/i^. We emphasize that our model parameters A and T, are adjusted to n(T, i/ij) at T = 1.1 
and have proven above to describe at the same time n(T,ifii) at T = 1.5, 2.5, 3.5 Tc. In so 
far, the agreement of our characteristic curve emerging at Tc with the transition line in [16] 
is quite satisfying. 

In Fig. 5 we also show the first two Roberge- Weiss transition lines (fat dashed curves 
characterized by [15] fi%/T^ = —T'^/T'^Tt^{2k + 1)^ for k = 1,2) at which thermodynamic 
quantities exhibit an analytic behavior at small temperatures, while the Roberge- Weiss 
transition represents a first-order phase transition (fat solid section) at sufficiently large T. 
In addition, the first center symmetry line is shown (dotted curve characterized by [15] 
IJ>%/T^ = —T'^/T^7!''^{2kY for k = 1). The repeated copies of these sectors for > 2 are not 
displayed in Fig. 5; they reside in the left bottom edge. 

Numerically, we find that the characteristic curve emerging at T = and the first 
Roberge- Weiss transition line cross each other at /Tc = 1.112 and {fi^Y/T^ = —12.214, 
whereas the lattice QCD simulations [16, 19] report T^/T^ = 1.095 and {fi^ f/T^ = -11.834. 
These tiny differences can hardly be resolved on the scale displayed in Fig. 5. For larger 
negative fi^ the characteristic curve is mirrored at the Roberge- Weiss transition line (see 
the section in the left top edge below the first Roberge- Weiss transition line in the left panel 



In the right panel of Fig. 5, we exhibit the phase diagram with the same notions as in the 
left panel, but with a linear abscissa ^b/Tc-, negative values of Hb/Tc are to be assigned to 



of Fig. 5). 
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purely imaginary chemical potential, while positive values correspond to purely real values 
ol jiB- As the coefficients in Peshier's flow equation Eq. (5) for real chemical potential obey 
^ for ^ and ^ for T ^ 0, the characteristic curves, including the one 
crossing the T axis at T^, approach the T = and Hb = axes perpendicularly. Deviations 
between QPM results for real hb and the polynomial fit become visible for /x^ > 330 McV. 
In addition, we exhibit the solution of Peshier's flow equation starting at Tc for a different set 
of QPM parameters (cf. section IV A) by the dashed line. For imaginary chemical potential, 
both results are indistinguishable, whereas for real /i deviations become visible for /ib/Tc > 2 
signalling again that small deviations in the imaginary chemical potential sector result in 
larger deviations in the sector of real fi. This makes predictions about the onset of possible 
deconfinement effects at small T and real /x difficult. 

C. Quark mass dependence 

In this subsection, we study the influence of a different approximation of the quark 
dispersion relation on the QPM results. In particular, we concentrate on the characteristic 
curve emerging at T = Tc. Apart from Eq. (3), we can approximate the quark dispersion 
relation by [28] 

co^ = k'^ + rnl + 2mqM+ + 2Ml (9) 

as for instance employed in [11, 29] with from Eq. (4). Changing the approximation of 
the dispersion relation cug demands a readjustment of the parameters in the effective coupling 
C^lT, /i = 0) in Eq. (6) in order to appropriately describe the lattice QCD data oin/T^ and 
causes changes in Peshier's flow equation (see Appendix A). The QPM parameters for using 
Eq. (9) adjusted to the continuum estimate of the n/T'^ data at T = 1.1 Tc read Tg — 0.976 Tc 
and A = 95 implying that the divergence in G^{T,iJ, = 0) is located at T = 0.987 Tc. With 
this new parametrization, the agreement between QPM and continuum extrapolated lattice 
QCD data is nearly as perfect as observed in Fig. 1 (at most 3% deviations). However, it 
indicates that Eq. (3) might be somewhat more suitable than Eq. (9) as quark dispersion 
relation. The influence on the characteristic curve emerging at T = Tc when employing 
Eq. (9) instead of Eq. (3) is negligible for imaginary chemical potential. For real i^b/Tc, the 
difference between both parametrizations is also very tiny and approximately 1.5% at small 
temperatures. 
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In addition, we can discuss the quark mass dependence of the found results by performing 
a naive chiral extrapolation niq — > 0. For imaginary chemical potential, quark mass effects 
turn out to be negligible independent of the specific quark dispersion relation used. For real 
Hb/Tc, quark mass effects are also small and visible only for very small temperatures. The 
differences between using = 0.2 T and = are less than 1% when employing Eq. (3) 
and at most 3% when employing Eq. (9). In both cases, decreasing quark masses imply a 
larger curvature of the characteristic curves and thus a smaller critical chemical potential at 
T — 0. A similar minor quark mass dependence with the same trend when decreasing rrig 
was found in lattice QCD simulations [30]. 

Another sensitive measure for the quark mass dependence of the net baryon density Ub 
is the chemical potential dependence of the chiral condensate < ipip > which arc related to 
each other by a Maxwell relation [19, 21]. Within the QPM, the /i dependence of < ipip > 
is fairly well described and will be reported elsewhere. We note that simply putting nig — 
but keeping the parametrization of ^^(T, = 0) fixed modifies n/T^ by less than 1% for 
the considered range of temperatures and chemical potentials. In principle, however, a 
general quark mass dependence of the parameters Tg and A in the effective coupling would 
be conceivable. Due to the minor effects observed, we restrict our further considerations to 
quark dispersion relation Eq. (3) in the following. 

D. fj, dependence of the quasiparticle masses 

In [18], the lattice QCD data [31] have been discussed with the goal to extract the relevant 
excitation modes from thermodynamic bulk quantities. The explicit fi dependence of the 
quasiparticle masses has been named BKS effect. In order to test the importance of the 
BKS effect on the found results, we omitted the /if/vr^ terms in the quasiparticle dispersion 
relations or flipped their signs though leaving the dependence of on /li unchanged. In 
fact, neglecting simply the term /if/n'^ in Eq. (4) (or Eq. (18) below), n/T^ is only affected 
for large //j, where the attenuation of by I/tt^ becomes smaller and the term proportional 
to iJ,f cannot be neglected compared to the term proportional to T^. This imphes that for 
larger T significant effects can only be seen at sufficiently large values of yUj. Note, however, 
that iii is restricted by /Xj < |T. Similar effects can be observed when flipping the signs 
in the asymptotic mass expressions. Nonetheless, thermodynamic self-consistency of the 
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QPM requires in both considered cases changes in Peshier's flow equation (5) rendering the 
coefficients b, and a^. according to Maxwell's relation. In such a thermodynamically 
self-consistent approach, we find our results to be indistinguishable from the QPM results 
exhibited in Fig. 1, i. e. found results seem to be rather independent of the explicit form of 
the /Xj dependence in the asymptotic mass expressions. However, a general dependence of 
the asymptotic masses on the chemical potential seems to be important. When neglecting /ij 
completely in the quasiparticle dispersion relations, thermodynamic self-consistency dictates 
also an independence of T in (and rrioo in Appendix A) which significantly changes the 
results: even though the almost linear behavior of n/T^ for small /ij can be reproduced, 
the pronounced curvature for T = 1.1 Tc at larger /ij cannot be obtained under such an 
assumption. The /ij dependence of n/T^ is further discussed in Appendix B. 

E. Scaling properties 

In [25, 32], a scaling of the ratio Ap/Ap^^ of the excess pressure in fiB direction was 
reported. Here, we find a similar scaling for the ratio ns/n^^ as depicted in Fig. 6, where 
denotes the Stefan-Boltzmann expression of the net baryon density. When considering 




1 1.2 1.4 1.6 1.8 2 1 1.1 1.2 

T/Tc T/Tc 



FIG. 6: Left: ratio n^/nl^ function of T/Tc for different imaginary and real baryo-chemical 
potentials. Dashed curves represent results for imaginary baryo-chemical potential, with |/xb/Tc| = 
0.6,1.2,1.8 from bottom to top, while solid curves depict corresponding results for real yU^, with 
l/f^s/Tcl = 0.6, 1.2, 1.8 in inverted order, i. e. from top to bottom. Right: zoom into the region close 
to Tc. 
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ns/n^ either for real or for imaginary chemical potential, in both cases, effects become 
visible only in the vicinity of Tc- In fact, for the baryo-chemical potentials considered in 
Fig. 6, the ratio UB/n^^ is found to be independent of //^ for T > 1.2 Tc. Furthermore, 
Ub/ti'^^ = 1 is approached only asymptotically, signalling the expected strong deviations 
from the free field behavior. Apart from the observed differences in the ratio between real and 
imaginary chemical potentials close to Tc, we find an interesting pattern: ub/ti^^ decreases 
with increasing baryo-chemical potential for real /ib, while the ratio increases in the case of 
imaginary chemical potential. This is partly caused by differences in ub/T^ between real and 
imaginary chemical potential which become smaller for increasing temperature (cf . Fig. 3) . 
But it is also related to different signs in for real or imaginary chemical potential. 

To be specific, in the case of real chemical potential, ns can be expanded into a Taylor 
series in powers ol jiB with expansion coefficients [31] Ck{T) — ^ ^ %Qj,}t)'^ ^ ■ "^^^ 
Stefan-Boltzmann expression for the net baryon density reads 

Even though this expression for nf^ is correct only for a massless ideal gas, while ub entering 

the ratio is evaluated for nig = 0.2 T, quark mass effects can safely be neglected (as discussed 
in section IIIC). The ratio UB/n^ reads 



riB ^ ^ ^ 



n 



SB 



Nf Nf \3TJ 



2C4 ^ 



0{^^%) . (11) 



In the hmit /j^b 0, the ratio approaches 2c2/Nf — ^x{T,iJ, — 0)/T'^ for Nf — 4. For 
small //b/(3T), i. e. for small /xb or large T, jiB effects become small, thus explaining the 
observed scaling. Furthermore, as 2c4 — C2/7r^ > for all temperatures T >Tc and remains 
approximately constant for T > 1.2Tc (cf. section IVB), a fixed ratio ns/n^^ requires 
increasing temperatures T for increasing ij,b, explaining the observed ordering in Fig. 6. 
Close to Tc, deviations between exact results and the Taylor series expansion of become 
larger with increasing /ib, such that the arguments presented here do not apply. 

In the case of imaginary chemical potential, ^^(T, i/ii) from Eq. (1) can be evaluated for 
small fii by expanding the trigonometric functions in powers of fJ,i/T yielding also a Taylor 
series expansion similar to the one in the sector of real chemical potential. Within this 
approach, we find 

nB{T, ijii) = i 0C2/Xir^ - ^C4/x- + ...^ , (12) 
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where 

1 d\p{T,iiii)/T' 



Ck{T) 



= Ck{T) . (13) 



m d{iii/TY 

Note that both, Ck and Cfe are real and Cfe, = for odd k. The Stefan-Boltzmann result 
for imaginary chemical potential reads 

n|^(r,z^.)=z(^/..r^-0^,^) (14) 

and the ratio follows as ^ — (fr)^ [2c4 — + 0{ii%). Similar to the considera- 
tions for real chemical potential, we observe a scaling with ^b/{^T) and in the limit //s — > 0, 
nB/n^ — > 2c2/Nf — \x{T,iJi — 0)/T^ for Nf — 4. For imaginary chemical potential, how- 
ever, the sign of the term proportional to ^\ is flipped, explaining the different ordering 
observed in Fig. 6, i. e. at fixed T, UB/n^ becomes larger with increasing hb- 



IV. COMPARISON WITH LATTICE QCD DATA AT /x = 



A. Pressure 



Via the QPM, we have access to both, real and imaginary chemical potentials. Thus, 
we can compare our results based on the lattice QCD data of [19, 21] with other lattice 
QCD calculations. In [33], a similar lattice setup for calculating the pressure at // = for 
Nf = A degenerate quark flavors with = 0.2T on a lattice with Nt- = A and = 16 
was considered, though employing an improved lattice action. These lattice QCD data [33] 
require also a proper continuum extrapolation. We apply a similar strategy as in section III A 
but now for the pressure, because its Stefan-Boltzmann limit is given in [33] for Nt — 4, and 
find — 0.839 as continuum extrapolation factor. The difference between and di^t in 
section III is maybe a consequence of the different lattice actions used in the simulations [33] 
and [19, 21] resulting in different cut-off effects on the data. 

In Fig. 7, the continuum estimated lattice QCD data [33] (squares) for p/T'^ as a function 
of T/Tc at // = are compared with the QPM using the parameters A and T, in G'^(T, fi = 
0) from section III A adjusted to n/T^ for imaginary chemical potential. The integration 
constant -B(Tc) adjusting the QPM value of p/T^ a,t /i — and T = Tc to lattice QCD reads 
B{Tc) = 2.567;^ using again = 163 MeV. As evident from Fig. 7, the general trend of 
p/T^ and the behavior at large T is reproduced. Nevertheless, p/T^ shows deviations of up 
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FIG. 7: Comparison of the QPM (solid curve, using parameters as in section III A) for the scaled 
pressure p/T^ as a function of T/Tc at = with the continuum estimate of the lattice QCD 
data [33] (squares). 

to 20% in the intermediate temperature region T ~ 1.5 Tc. The same deviation pattern was 
already discussed in [14, 25, 29]. In fact, it seems to be a general feature, that fits to lattice 
QCD data in the sector of zero (non-zero) chemical potential underestimate (overestimate) 
the according results in the sector of non-zero (zero) chemical potential. 

Considering, instead, an independent adjustment of the QPM parameters to p/T^ at 
= 0, the comparison of the QPM with the continuum estimate of the lattice QCD data is 
exhibited in Fig. 8 (left panel). We find an impressive agreement when adjusting = 0.91 
A = 16 and B{T^ — 1.25 with — 163 MeV. With this new parametrization, we evaluate 
the net quark number density for imaginary chemical potential. The results (solid curves) 
are shown in the right panel of Fig. 8. At constant fJ^i/Tc, we find increasing deviations for 
decreasing temperatures, in particular close to the Roberge- Weiss transition, even though 
the change in slope close to /ic is qualitatively still reproduced. Furthermore, using this 
set of parameters results in the dashed curve in the right panel of Fig. 5 as characteristic 
curve emerging at T = Tc for real hb- With respect to the observed deviations at T = 1.1 Tc 
between the lattice QCD data from [19, 21] for n/T^ and the QPM with parameters adjusted 
to p/T'^ (see right panel of Fig. 8), it is surprising that this characteristic curve agrees so 
well with the one considered in section III B. 



17 



B. Taylor expansion coefficients 

Having the QPM parametrizations employed in Figs. 1 and 8 at hand, we can discuss 
their influence on the Taylor series expansion coefficients Cfc(T) defined in section HIE for 
Nf = 4 similar to studies for Nj = 2 in [14]. In Fig. 9, we exhibit C2(T) and C4(T): 
C2(T) = ^x{T,fi = 0)/T^ shows some deviations between both parametrizations whereas 
C4(T) agrees fairly well with visible deviations only in the very vicinity of T^. In addition, 
we observe that smaller co(T) = p{T,fi = 0)/T^ (cf. Fig. 7 and the left panel of Fig. 8) 
implies smaller C2(T) as already pointed out in [29]. As already mentioned in section III A, 
n/T^ depicted as a function of /ij/T shows almost no temperature dependence for T > 1.5 T^. 
This is mainly due to the fact that C2(T) (upper solid curve in Fig. 9) exhibits also a rather 
negligible temperature dependence for larger T. Furthermore, C4(T) is sizeable only close 
to Tc and approaches its Stefan-Boltzmann limit l/vr^ for T > 1.2 Tc. Considering ub/T^ 
in terms of a Taylor series expansion up to order the results for real and imaginary 

chemical potential differ only in the sign of the cubic term which is oc C4/i^. Thus, the net 
baryon density evaluated for real or imaginary chemical potential deviates only for larger 
chemical potentials and close to as evident from Fig. 3. Cq{T) (not exhibited) deviates 




FIG. 8: Left: comparison of the QPM (solid curve) for p/T^ as a function of T/T^ at /i = 
employing readjusted parameters (see text) with the continuum extrapolated lattice QCD data 
(squares) as exhibited in Fig. 7. Right: comparison of the QPM (solid curves, using the readjusted 
parameters) for n/T^ as a function of /ij/Tc for different T with the continuum extrapolated lattice 
QCD data (symbols) as exhibited in Fig. 1. 
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FIG. 9: Taylor series expansion coefficients C2(T) (upper curves) and C4(T) (lower curves) as 
function of T/T^ for Nf = 4 employing the different parametrizations from Fig. 1 (solid curves) 
and Fig. 8 (dashed curves). 

significantly from zero only for temperatures very close to but can become of the same 
order of magnitude as c^iT) at T = Tc. 

Finally, we mention that the parametrization in section IV A, optimized for reproducing 
the Nf = 4 lattice QCD data of [33], can also be used to describe the Taylor coefficients 




FIG. 10: Left: comparison of lattice QCD data [31] for Nf = 2 for C2(T) (squares) and C4(r) 
(circles) as a function of T/Tc with the QPM for Nf = 2. Solid curves represent results applying 
QPM parameters as in Fig. 1 adjusted to n/T^ while dashed curves represent results applying the 
parametrization of Fig. 8 adjusted to p/T^ at /x = 0. Right: ratio C4/C2 as a function of T/Tc for 
both parametrizations. 
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C2,4:{T) from lattice QCD [31] for Nf — 2 as exhibited in the left panel of Fig. 10 (dashed 
curves) . In doing so, we keep A and Tg adjusted to Nf — A lattice QCD data fixed and change 
merely from Nf — A to Nf — 2 in the thermodynamic expressions of the QPM. In fact, we 
find deviations of about 10% close to Tc and less than 5% for T > 1.2 Tc between QPM 
and [31] for C2(T). While this coincidence might be accidental, one could also argue that 
the quasiparticle model catches correctly the flavor dependence. In contrast, employing the 
parametrization from section III A gives a pattern resembling Fig. 9. Even though deviations 
between both parametrizations are obvious, the ratio C4/C2 is rather insensitive with respect 
to the employed parametrization for T > 1.2 Tc approaching l/(27r^) as shown in the right 
panel of Fig. 10. 



V. SUMMARY AND DISCUSSION 



In summary we extend our effective quasiparticle model and compare it with lattice QCD 
data for purely imaginary chemical potential. Despite the fact that our phenomenological 
model does not exhibit the Roberge- Weiss periodicity of full QCD, it is able to describe the 
available lattice QCD data [16, 19, 21] impressively well. In particular, the drastic change 
in slope of n/T^ close to the critical chemical potential fJ^c/T = 7r/3 of the Roberge- Weiss 
transition can be described. This is entirely due to the BKS effect [18], i. e. a consequence 
of chemical potential dependent quasiparticle masses. A thermodynamically consistent in- 
vestigation of the importance of the //j dependence in the quasiparticles' asymptotic masses 
shows that the found results are independent of the chosen explicit form of the /ij depen- 
dence. Nonetheless, the pronounced structures cannot be reproduced when the quasiparticle 
masses would be completely independent of /ij. In this respect, the fii dependence imple- 
mented in the model is confirmed. Another evidence is the comparison of the QPM result 
for the characteristic curve emerging at T = Tc with the phase transition line evaluated in 
lattice QCD simulations [16]. For the Roberge- Weiss transition, we find critical values of 
temperature and baryo-chemical potential close to the ones given in [16]. The successful 
comparison points to the correctness of Peshier's flow equation as a tool for transporting 
information from /x = to non-zero /i which is of particular importance for the knowledge 
of the equation of state at larger baryon densities relevant for CERN-SPS and upcoming 
FAIR. 
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With the found QPM parametrization describing lattice QCD data in the sector of purely 
imaginary chemical potential at hand, we can also compare with an independent set of 
lattice QCD data [33] obtained at = 0. We find some deviations for the pressure in the 
intermediate temperature region which might account for the different lattice actions used 
in the calculations but could also signal to some extend a disagreement of results obtained 
at /J, — and non-zero /i, as already discussed in [14, 29]. In this context we emphasize 
that the comparison of thermodynamic models with lattice QCD data is hampered by the 
lacking systematic continuum extrapolation of the latter. 

Ab initio it is not clear whether the assumed quasiparticle excitations represent the proper 
description of QCD thermodynamics also in the region close to T^. The success of the present 
comparison lends some credibility into the picture of quasiquark excitations with a mass gap. 
This in in line with findings in [34], where also a striking deviation from the pertubative 
excitation pattern close to Tc has been found. Nevertheless, it would be premature to claim 
that the strongly coupled hot quark-gluon medium is entirely described by the presently used 
quasiparticles. For instance, excitation modes like plasminos and longitudinal gluons are not 
included in the model. Furthermore, finite width effects of the quasiparticles and Landau 
damping are neglected. One should keep in mind that thermodynamic bulk properties are 
sensitive essentially to excitations with hard momenta, i.e. /c ~ T, /x. There may be a variety 
of soft and ultra-hard excitations rendering the picture of the strongly coupled quark gluon 
plasma into a much more involved scenario, in line with the complexity of QCD. 

Our model is far from being an ab initio calculation as attempted in [35]. But, in particu- 
lar, the fiexibihty of the introduced effective coupling allows for curing possible deficits in 
the dynamical degrees of freedom. Apart from that, the model is highly non-perturbative as 
it can be formulated in terms of an infinite series of powers in the coupling, though, making 
contact with perturbation theory, as the first terms coincide with perturbative QCD and 
asymptotically approaches the running QCD coupling. 

Finally, we remind the reader that we consider here a fairly special case of four degenerate 
quark flavors. Despite of the known sensitivity of particular features of QCD on the flavor 
content, some scaling properties of thermodynamic bulk properties may be useful for an 
orientation in thermodynamic state space. 
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Appendix A - Flow equation for imaginary chemical potential 

The QPM pressure for imaginary chemical potential reads 

p{T, ini) = Pa{T, i^i) - B{T, ijii) (15) 

with partial pressures 

Pq{T, ijii) = / dkk^ (In [1 + e^*'^'-"^)/^] + In [l + e(-*^>— ^)/^] ) , (16) 

Pg{T, If,,) = -^T / dke In [1 - e-'^^/^ (17) 
Jo 

for quarks and gluons, respectively, where dq — 2NcNf and dg — — 1. The quasiparticle 
dispersion relations read cu^ — k^ + rUg + 2M^ = k^ + with given in Eq. (4) and 
cUg — k^ + ml with asymptotic mass 

ml = ^ ([27V. + Nf]T' - ^Nff,^ G'{T, z/.,) . (18) 

Assuming that all T and /ij dependence of the function B is encoded in the asymptotic 
mass expressions Mqo and moo, thermodynamic consistency is fulfilled from the stationarity 
conditions dp/dM^ = dB/dM^ and dp/dm1 = dB/dm1 [36] such that entropy density s 
and net quark number density n are obtained from standard thermodynamic relations. The 
purely real result for s — Sq-\- Sg reads 

Sq{T, = (^^^ [ft + /.1 - ^^^^ [ft " /.1 ) ' (19) 



(20) 

Jq ujg e'^al^ - 1 



where = (^f.{u^,^ii^i)/T ^ ^^-i g^^d n is given in Eq. (1). The quasi- linear partial differential 
equation Eq. (5) to be solved for G^(T, i/Xj) follows from Maxwell's relation 

ds d^p d^p dn 



d{ifMi) d{iHi)dT dTd{tfMi) dT ' 
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(21) 



where the exphcit derivative terms cancel each other leaving 



dn dM^ _ dsq SM^ ^ dsg drri^ 



(22) 



dMl dT dMld{iiii) dm1d{iiii) 

Omitting the overall factor of i in dn/dM^, dM^/d{iiii) and dm'^/d{iiJ,i), the coefficients 
of Eq. (5) read 

Cf 



^TG^ + 2mga ] h 



27r2 



C 



f 



rp2 



II, 



12 



l[2N, + Nf]T'- 



NrN: 



(23) 
(24) 
(25) 



where Cf — {N^ — 1)/{2N^ and the integral expressions exphcitly read 



Ii 
I2 
I3 



dn 



27r2T 

dq 
27r2T 



dk 



P (e"«/^sin(^,/T) - e3-«/^sin(/i,/T)) 



Uq (e2^./r + 2e^^/T cos{fii/T) + 1)^ 



00 T 2 

dk- 



1 



'4^.2 
^3 



A;2 + m^) 



1 + 



dk 



A;2/ 2e^^/^cos(/ii/T) + 2 
^ ye^^^lT" + 2e^-il^ cos{fii/T) + 1 



T (e^s/^ - 1) 



(26) 
(27) 



+ Ml) {26^^"/^ cos{fii/T) + 4e2-«/^ + 2e"^/^ cos(/Xi/T)) 



'4^,2 
^3 



2CUgT 



(gS-'./T + 2e-«/^cos(/ii/T) + 1)' 



/x^ (e^-'/^sin(^^/T) - e^^"/^ sm{fi,/T)) \ 
T {e^^./T ^ 2e^^/T cos{iii/T) + if ) 



(28) 



The term in b proportional to a stems from assuming temperature dependent quark masses 
niq — aT as employed in some lattice QCD performances, e. g. [19, 21, 31, 33, 37]. 

The quark number susceptibility x Eq. (7) is found to be symmetric when replacing 
/Ii by —Hi because the same holds true for dG'^/dfii. From Eq. (5) we find 

b ar dG^ 



a„. dT 



(29) 



For the individual expressions entering Eq. (29) we find Ii — > — Ii, I2 — > I2 and I3 — > I3 
for /Ii — > —/Ii such that b —b, ax —or and a^. — > a^^.. In addition, a Taylor series 
expansion of G^{T,ifj,i) in powers of /Xj consists only of even powers in /Xj [13, 14] such that 
dG^/dT is symmetric under /ij — > — /Xj. In the limit /li — > 0, we find dG^/d/ii — > as 6 — > 0, 
qt ^ but a^. and dG^ /dT remain non-zero. 
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When employing the quark dispersion relation Eq. (9), the coefficients of the flow equation 
render to 

/ 



c 



f 



+ ma 



Cf 



^ 8 (r2 - 



2TG^ + 2mqa + 2aM+ 



Ii 
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—r + ruq 
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Appendix B - Parametrizing the Hi dependence 

In section III, we found a general dependence of the quasiparticle dispersion relations on 
temperature and chemical potential to be of utmost importance for the successful description 
of lattice QCD data. This shall be illustrated in some more detail by considering the net 
quark number density in Eq. (1) of an ideal gas with dispersion relation uj'^ — k"^ -\- M^. 
In principle, thermodynamic self-consistency demands cither a dependence of M on both, 
T and /Xj, or neither a T nor a jii dependence of M . In the latter case of constant M, 
we adjust M — 0.21 GeV in order to describe the continuum extrapolated lattice QCD 
data of n/T^ (cf. Fig. 1) at T = 1.1 Tc for small Hi/T^. The corresponding QPM results for 
T — 1.1, 1.5, 2.5, 3.5 Tc are then exhibited in the left panel of Fig. 11 (dashed curves). At T = 
1.1 Tc, we find increasing deviations from the lattice QCD data for jii/Tc > 0.66, in particular 
in the vicinity of the Roberge- Weiss critical chemical potential fi^, where the pronounced 
curvature cannot be reproduced. This was already discussed in [21] by considering the ratio 
n{iii) / n{iii) free Signalling clear deviations of the lattice QCD data from a free (ideal) gas 
behavior. Furthermore, by increasing T but keeping M fixed, the description of the lattice 
QCD data becomes less and less accurate for smaller Hi/Tc suggesting a general dependence 
of M on T. Readjusting M individually for each temperature, ignoring for the moment being 
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FIG. 11: Left: comparison of ideal gas results for n/T^ as a function of /Uj/Tc employing either a con- 
stant mass parameter M = 0.21 GeV (dashed curves) or readjusting M/T = 1.17,0.90,0.81,0.77 
for T = 1.1, 1.5, 2.5, 3.5 Tc (soHd curves from top to bottom) with the continuum extrapolated 
lattice QCD data (symbols) as exhibited in Fig. 1. Right: comparison of found M/T (squares) 
as a function of T/Tc with the asymptotic quark mass Mqo/T of the QPM at /x = employing 
Ts = 0.96 Tc and A = 56 as in section III A. 

thermodynamic self-consistency, the results are depicted by solid curves in the left panel of 
Fig. 11. The found scaled mass parameters M/T for the temperatures considered here are 
exhibited in the right panel of Fig. 11 (squares) and compared with the scaled asymptotic 
quark mass M^o/T of the QPM at /i = (solid curve) employing the parametrization of 
section III A. Both results agree fairly well, indicating that non-zero chemical potential 
effects are tiny for small /ij/T but become sizeable close to fJ^dT) as also visualized in 
Fig. 12. In Fig. 12, the scaled asymptotic quark mass M^/T of the QPM is exhibited as a 
function of fii/Tc for constant T using the QPM parametrization of section III A perfectly 
describing n/T^ in Fig. 1. For increasing T, M^/T shows decreasing sensitivity on /i, while 
non-zero chemical potential effects become important close to /Xc = fT"- 



25 




FIG. 12: Asymptotic quark mass M^olT of the QPM (solid curves) as a function of fii/Tc for 
T = 1.1, 1.5, 2.5, 3.5 Tc (from top to bottom) employing the QPM parameters from section III A. 
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